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Abstract 

We show that the different ways of deriving the Heavy Quark Effective Theory 
(HQET) lead to equivalent theories. The equivalence can be established through a careful 
redefinition of the field variables. We demonstrate the equivalence to order ^ in the 
presence of a constant electric field. 
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I. Introduction: 

In recent years, there has been a lot of work in the context of Heavy Quark Effective 
Theories (HQET) [1-6]. When a quark is heavy, it can be effectively described by a two 
component spinor and there are essentially two different ways of obtaining an effective 
theory for such a quark. The most logical way to take the heavy quark mass limit, in 
some sense, is through the Foldy-Wouthuysen transformations [7,8] where one essentially 
diagonalizes the quark Lagrangian such that the upper and lower two component spinors 
decouple after which one restricts to the positive energy two component upper spinors. 
(In the functional language, one would simply integrate out the lower component spinors 
after diagonalization and absorb the result into the normalization factor.) The other 
approach to an effective theory is the traditional way [4,9] where one decomposes the 
four component spinor into a "large" and a "small" two component spinor. One then 
eliminates the "small" component through its equation of motion to obtain an effective 
theory for the "large" component spinors. (In the functional language, one would integrate 
out the "small" components and absorb the result into the normalization factor of the path 
integral.) While the two methods give the same identical theory in the lowest order in ^, 
their structures appear to be different at higher orders. This disagreement has justifiably 
raised concern in the recent literature [10] mainly because higher order corrections to 
various physical processes are currently being calculated [11] using the theory resulting 
from the traditional method of eliminating the "small" components. 

The discrepancy between the resulting theories in the two approaches is nothing new. 
It was already noted in connection with an electron interacting with an external electric 
field [12] as well as in the context of the Tamm-Dancoff method in nuclear physics [13], 
that eliminating the "small" component naively leads to a nonhermitian Hamiltonian. 
(In the case of the electron interacting with an external electric field, the lowest order 
manifestation of the nonhermiticity is in an imaginary electric dipole moment.) To the 
best of our knowledge, the resolution of this puzzle was first proposed in ref. 14 where 
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it was noted that a renormalization (redefinition) of the large components is essential for 
the hermiticity of the Hamiltonian. In the case of the Dirac electron interacting with an 
electric field, the equivalence of the two approaches was thereby demonstrated to order 



It is in general believed that an appropriate field redefinition will lead to an equiva- 
lence of the two approaches. In this note, we wish to demonstrate up to order ^ that 
both the methods indeed give the same theory with appropriate field redefinitions. Both 
the methods have their advantages and disadvantages and we comment on this in the 
conclusion. The equivalence of the two approaches is quite important and so is the un- 
derstanding of the field redefinition since otherwise the higher order corrections calculated 
with the naive effective theory may not represent the true physical effects. In sec. II, we 
review the known results up to order ^ and try to bring out the necessity for a field 
redefinition. In sec. Ill, we show that, in the case of the free theory, the two approaches 
lead to the same effective theory upon field redefinition. We then demonstrate to order 
^ that in the presence of a constant electric field, the two approaches also give the same 
effective theory through a careful redefinition of fields. In sec. IV, we show the equivalence 
of the two approaches in the functional approach and present our conclusions in sec. V. 

II. The Problem of Imaginary Dipole Moment: 

In this section, we consider a Dirac electron interacting with a constant, external 
electric field. While the discussion can be carried out equally well in the first or second 
quantized language we will follow a first quantized approach for simplicity and clarity. Our 
discussion in this section will follow closely the work in ref. 14. The Dirac equation in the 



[14]. 



present case has the form (we use Bjorken-Drell metric, 7' 
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(2.1) 



with 



— * — * 

E = VAq = constant 



(2.2) 
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If we let 

^ ^ e-^"^V (2.3) 

then Eq. (2.1) takes the form 

(^7° (^0 + ieAo) + if • V - m (l - 7°)) V = (2.4) 

Introducing the "large" and the "small" components as (In the current terminology of the 
subject our entire discussion will be with the choice v^^ = (1, 0, 0, 0) for simplicity.) 

-(::) 

we note that the Dirac equation (2.4) separates into two equations 

= {ido - eAo + 2m)"' (^-ia • v) V'l = (2-6) 

{ido - eAo) iJL = -iB ■ VtPs (2.7) 

Upon substituting Eq. (2.6) into Eq. (2.7) and expanding up to order the equation 
for the "large'' component takes the form 

(»„ - eA„) = (- J- - . X V) + E . V + O (ij)) (2.8) 

It is the last term on the right hand side in Eq. (2.8) which represents an imaginary 
electric dipole moment and arises naturally as a consequence of eliminating the "small" 
component spinors. While this may be puzzling, it's origin is not hard to understand. 
As explained beautifully in ref. 14, one can view the process of eliminating the "small" 
component also equivalently as finding a transformation which will take 

'■=(:;) 

Such a transformation is generated by the matrix 



It is clear that this transformation is not unitary and as a result the Hamiltonian 
does not remain hermitian under such a transformation. Another manifestation of the 
transformation not being unitary is to note that under such a transformation, the norm of 
the state is not invariant. In fact, let us note that 



(2.10) 



Therefore, if 

Jd^xtp^ip = l (2.11) 

then. 



J d^x ipli^L 1 (2.12) 



It is clear, however, that the norm can be maintained (state will be normalized) if we 
redefine 

V^z.= (1 + A^A)'^'V'L (2.13) 

As is described in ref. 14, this redefinition also restores hermiticity of the Hamiltonian. 
(The norm will, of course, be time independent.) In fact, let us note from Eqs. (2.6) and 
(2.13) that if we define up to order 



1 



or, V'L ^ ( 1 + 77^ \ 



(2.14) 



8m2 

then Eq. (2.8) would take the form 

_ e^„) = (- J- - . X V) + O (ij)) (2.15) 

The imaginary electric dipole term has now been absorbed into the field redefinition and 
as we will see in the next section this is exactly the same equation which we would get 
through the Foldy-Wouthuysen transformation up to order 
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III. Equivalence up to Order 

As a warm up, let us consider a free Dirac particle satisfying 

(i7°9o + Z7 • V - m) V = (3.1) 

The idea behind the Foldy-Wouthuysen transformations [7] is to block diagonalize the 
Dirac operator such that the upper and the lower two component spinors are decoupled. 
It is well known [7] that under the unitary transformation 



where 



tanh-MM) (3.3) 



m I |v| \ "7, 

the Dirac operator transforms into the diagonal form 

(z7°ao + Z7 • V - m) ^ U{i^^do + i^-V-m)U^ 

= i7%-(m2-V2)'/' 

Thus the transformed equation has the form 



Writing 

( 

V'2 



(3.4) 



(^ij^'do - (m^ - V^) '^'^ V' = (3.5) 

If we further let 

then Eq. (3.5) will take the form 

(^z7°ao + m7° - (m^ - V^) '^^^ V' = (3.7) 



^' = I (3.8) 
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and restricting to the upper two component spinor, we obtain from Eq. (3.7) 

(^ido + m- (m^- V^) ^^^^ Vi = (3.9) 

A power series expansion of Eq. (3.9) in ^ would then give the effective dynamical equation 
in the Foldy-Wouthuysen approach. 

In contrast, let us note from Eqs. (2.6) and (2.7) that in the absence of interactions 

il;s = A^Pl = {ido + 2m)-^{-i& ■ V)V'l (3.10) 
idQipL = -io ■ Vips = ■ ^{ido + 2m)" V • VipL = -{ido + 2m)~^V'^ipL (3.11) 

The equation for the "large" component, Eq. (3.11), of course, does not at all resemble 
Eq. (3.9) obtained through the other approach. Let us also note that if we define 

1 /2 

V^L= (1 + AU)'^'Vl= (l-(iao + 2m)-'v') (3.12) 

as we should to maintain normalization, the form of the equation does not change 

ido^L^ -{idQ + 2m)~^V^^L (3.13) 

This is mainly because the redefinition factor in Eq. (3.12) is field independent, commutes 
with as well as V and that Eq. (3.11) is linear in ■0/,. However, we note that Eq. (3.13) 
can also be written as 

{ido + {ido + 2m)-^V2)V^L = (3.14) 
Therefore, if we further redefine 

then Eq. (3.14) will take the form 

{^do + m-{m^ - V^) ^'^^ 'ipL ^ (3.16) 
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which is, of course, what we have in the Foldy-Wouthuysen approach (see Eq. (3.9)). 
The point of this exercise is to note that even after the spinor is normahzed properly, 
there stiU remains an arbitrariness up to unitary transformations [15] and this is crucial in 
establishing the equivalence of the two approaches. Alternately, we note that Eq. (3.14) 
can also be written as 

{ido + 2m) (^ido + m+ (m'^ - V^)^^^^ (^ido + m - (m'^ - V^)^^^^ V'l = (3.16') 

Consequently, if we redefine 

ido + m + (m'^ - W'^^ j i>L (3.16") 
then the dynamical equation would take the form 



{^do + m - (n? - V^j ^ ^ = 



which is the equation obtained through the Foldy-Wouthuysen approach. It would seem, 
therefore, that at the level of the equation of motion, there is an arbitrariness in the choice 
of field redefinition. Furthermore, the redefinition in Eq. (3.16") would seem to imply a 
change in the norm of the state. However, in the next section, where we will show the 
equivalence of the two approaches for the free theory in the functional approach, it will 
become clear that such a field redefinition is essential for the matching of the functional 
determinants. 

Let us next consider a Dirac particle interacting with a constant, external electric 
field. This is chosen to avoid the unnecessary complexities associated with a generalized 
non Abelian gauge potential which do not provide any new insight into the problem. The 
equation of motion is 

(i^^{do + ieAo) + i^-V -7x^11^ = ^ (3.17) 

with 

^ = VAo = constant (3.18) 



In the presence of interactions, the Foldy-Wouthuysen transformations have to be imple- 
mented order by order. We refer the readers to the hterature [16] for details and merely 
quote here the result that the equation satisfied by the upper two component spinor in the 
present case, up to order after a phase redefinition of the form in Eq. (3.6) is given by 

Note here that up to order it is the same equation as in Eq. (2.15). 

To analyze the effective equation in the other approach, let us recall from Eqs. (2.6) 
and (2.7) that in the present case 

ips = AipL = - {ido - eAo + 2m) a • VipL 

{ido — eAo) i/jl = —i^ ■ V^'s = ■ V {ido — e^o + ^m)'^ a ■ VipL 
The normalization factor, in this case, 

1 /2 

{l + A^A)^^^=(l-a-V {ido -eAo + 2m)~^a- v) (3.20) 

is field dependent and involves the Dirac operator. Expanding to order we note that 
the redefined field will be given by 

i^L = {l + A^A)'^^i;L 
or, ijL = {l + A^A)-'^^ipL 

3 ^ 2 - 1 - 3 - (3-21) 

^ V'^ &-V{ido-eAo)&-V 



128m5 

^ ff-V{ido-eAo)ff-VV^ + o{^-^ 



128m5 



Substituting this back into the equation for the "large" component we obtain after some 
algebra 



(ido - eAojipL 



■ 



2m 

e 

32m4 



Aw? ^ ' IGm^ 

(2^-V-3ia- (Ex V))V^- 



3 



8m3 



16m^ qo^5 V / 



32m5 



16m^ 

+ :i V2(iao - e^o) + ^ • V + i<j • X V)) (^^0 - eAo) 



8m2 ^ Sm^ 



3e 



1 



16m4 



+ 



(3.22) 



We note here that since the normalization factor involves the Dirac operator, Eq. 
(3.22) has to be solved iteratively. This is crucial because this implies that one must keep 
terms consistently up to any given order. We note that this is a new feature not present 
at order Iterating Eq. (3.22) up to order gives 



(i9o - eAojipL 



le 



— - 
2m 4m2 



1 



-(£'xV)-£3(VT + 



8m^ 

1 '^2n3 



+ 



32m5 



16m^ 



2^2 



(3.23) 



(^. v)Vo 



1 



rrr 



Comparing with Eq. (3.19), we see that the two equations almost agree - in fact, only 
the coefficients of the last two terms are different. However, as discussed earlier, we also 
recognize that i/^l is unique only up to a unitary transformation. Taking advantage of this, 
we let 



25e 



e 192m 



(3.24) 
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With this redefinition, we note that Eq. (3.23) takes the form 
[ido - eAo)V'L = 



2m Am? ^ ' 8m^ ^ ' 8m^ 



ff.(£xV)V^--Lj(V^)' + -5^£^V^ (3.25) 



IGm*^ 16m^ 48m^ 



24m5 

This is exactly the same as Eq. (3.19) which is what we would obtain in the Foldy- 
Wouthuysen approach. Thus, the two approaches give the same effective equation up to 
order ^ with appropriate field redefinition. 

IV. Functional Equivalence: 

In this section, we will show the equivalence of the two approaches in the functional 
method. Let us start with a free Dirac theory described by 

C = ^ {ip - m) %l) (4.1) 

If we make the phase transformation 

^ ^ e-^""V (4.1') 
then the Lagrangian will take the form 

£ = V'(«^-m(l-7°)) V (4.2) 
The generating functional in this case is given by 

Z = N J V^Vi/; e*/'^'"^ ^ (4.2) 
In the Foldy-Wouthuysen approach if we make a unitary redefinition of the field variables 



2 
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with 

jMU^tanh-fM) (4.4) 
\^ m y |V| \ / 

then the Lagrangian diagonahzes and takes the form 

C = ipUido + m-(m''- V^y^^Ui + 4(ido + m+(m^- ^^^^^^2 

(4.5) 

+ surface terms 

The functional measure does not change under the unitary redefinition of the fields in Eq. 
(4.3). Consequently, the generating functional has the form 

= iV[det {tdo + m+{m^- V^^^)] j p^tp^^gi/d^4^1(^^o+m-(m^-v^)^/^)^0 
= 7V[det {ido + m + (m^ - V^)^/^)] Zpw 

(4.6) 

We note here that the propagator in the Foldy-Wouthuysen approach would have the form 

E + m- (P + 777,2)1/2 

Clearly, this has a pole at 

E= (P+m2)i/2-m (4.8) 

with unit residue. 

On the other hand, if we write 

-(::) 

then the Lagrangian of Eq. (4.1') will have the form 

C = ijjiidoipL + i^iiff ■ Vips + ^li^ ■ VV'L + {ido + 2m) i;s (4.10) 

Since the Lagrangian is at the most quadratic in the fermions, we can integrate out the 
"small" component and obtain 

= N[det{ido + 2m)] J V^^^V^jl / '^'^^l(^«o+(i9o+2m)-iv2)^,, (4.11) 
= N[det{ido + 2m)] Zconv 

12 



From the structure of the Lagrangian in Eq. (4.11) we note that the conventional (naive) 
propagator for the "large" component will have the form 

^CONV = — 7^ (4-12) 
We note that the propagator has poles at 

1/2 



E 



= ±(^P + m2) -m (4.13) 



The location of the positive energy pole, in which we are interested, coincides with that 
of the other approach. But we note that the residue at the positive energy pole is 



_i_ f I ■^\^ 

,rn /o • This suggests that the field '0r needs to be redefined. Note that if we 
redefine (compare this with Eq. (3.16")) 

(1 /2\ ''"/^ 
z^o + m + (m^ - V^) j (4-14) 

then the generating functional will have the form 
Z = N[det{ido + 2m)] j PVl^V'Le'/^'^ ^l{id,Hid.+2m)-^v^)^^^ 
= N[det{ido + 2m)][det(iao + 2m)] 

X [detiido + m + (m^ - V^)!/^)] J v^[Vi^Le' I ^fl^id^+m-im'-v')'/')^, 

= N[det{ido + m + (m^ - V^/^)] Zconv 

(4.15) 

Comparing Eqs. (4.6) and (4.15), we conclude that 

Zpw = -^CONV (4-16) 

Namely, the two approaches lead to equivalent theories but only after appropriate field 
redefinition. For an interacting theory, Foldy-Wouthuysen transformations do not have a 
closed form. Therefore, the equivalence in the functional approach has to be shown order 
by order. An an example, let us consider the Dirac particle interacting with a constant 

13 



electric field up to order The Lagrangian with the phase redefinition of Eq. (4.1') is 
given by 



£ = ^ (^i/ + ieAo) + Z7 • V - m (1 - 7°)) V 



(4.17) 



In the Foldy-Wouthuysen approach, if we diagonalize the Lagrangian up to order it 
will have the form (FW transformations are unitary and hence the functional measure does 
not change.) 

C {^^\dQ + ieAo) - m(l - 7^) 



1 - 



2m Am? 



Writing if)' as in Eq. (4.3), the Lagrangian density has the form 



(4.18) 



£ = ( z^o - e^o + — + ^ a ■ (E X V) + O (-3 ) ) 



+ V'|(iao-e^o + 2m-^V^ + £^a-(^xV)+0(^))^2 



(4.19) 



The generating functional, therefore, has the form 



Z = N 
= N 



1 ie _ „/ 1 



det{ido-eAo + 2m- — W' + ^a-{ExV) + 0{-^) 



X 



J PV'JPV'ie*/'''" ^liido-eAo+^ v^+_^a=.(Exv)+o(^))Vi 



« ( 2m(l + ^ (>a„ - eA„) " 5^ + 0(^3)) 



4m^ 



(4.20) 



On the other hand, the Lagrangian in terms of "large" and "small" components has 
the form 

£ = (ido - eAo) ipL + ^[^a ■ Vips + i^U^ ■ V^'l + -^s (^^o - eAo + 2m) ips (4.21) 
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Integrating out the "small" components we have 



N 



det ( 2m(l + ^{ido - eAo)) 



X e 



TV 



1 



det 2m (1 + — (ido - eAo)) 



-^CONV 



If we now redefine (as in Eq. (2.14)) 



^ = ( 1 - ) 



then the Lagrangian in Eq. (4.22) would take the form 



/:coNV = ^1 ( - e^o + 2^ + £^ a ■ X V) + O ) ) V'L 



(4.22) 



(4.23) 



(4.24) 



Furthermore, with the Jacobian arising from the redefinition in Eq. (4.23), the generating 
functional would become 

-\ 2 



Z = N 



det 1 - 



1 

8m2 



det (^2m(l + ^ [ido - e^o))^ 



= N 



iet{2m(l + l-(ia„-eA„)-J-,V-+0(±)) 



^CONV 



Comparing Eqs. (4.20) and (4.25), we conclude that 



ZfW = ZcONY 



(4.26) 



This establishes equivalence of the two approaches up to order ^ in the presence of a 
constant electric field. Once again, the significance of the field redefinition cannot be over 
emphasized. It is tedious, but as the discussion of the earlier section shows, the equivalence 
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can be carried out to any given order with appropriate field redefinitions. The higher order 
equivalence, however, can only be shown through a careful iterative procedure. 

V. Conclusion: 

In this note we have shown up to order ^ that the two ways of obtaining an interact- 
ing heavy quark effective theory yield the same result. From a practical point of view, it 
is the Foldy-Woiithuysen method that is simpler for the derivation of the effective theory. 
However, since it involves field dependent unitary transformations, unless the fields vanish 
sufficiently rapidly, it is conceptually unclear whether the S-matrix elements will remain 
the same [17]. The traditional method of eliminating the "small" components, on the other 
hand, is quite tricky. Here the "large" component fields must be renormalized (redefined) 
and this must be carried out carefully in an iterative manner. This is essential since oth- 
erwise the calculations may not represent the true physical effects. This procedure of field 
redefinition is quite tedious but has the virtue that the existence of the S-matrix elements 
is never in question in this approach since the asymptotic form of the wave function re- 
mains unchanged in general. Thus, in some sense, one can view the two approaches to be 
complementary and equivalent. 
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